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Resume 

In this paper, we present a determinist Jordan normal form algorithms 
based on the Fadeev formula : 

(A-/ - A) ■ B{\) = P(A) ■ I 

where -B(A) is (A ■ / — A)'s comatrix and -P(A) is A's characteristic poly- 
nomial. This rational Jordan normal form algorithm differs from usual 
algorithms since it is not based on the Probenius/Smith normal form but 
rather on the idea already remarked in Gantmacher that the non-zero col- 
umn vectors of -B(Ao) are eigenvectors of A associated to Ao for any root 
Ao of the characteristical polynomial. The complexity of the algorithm 
is 0{n'^) field operations if we know the factorization of the character- 
istic polynomial (or 0{n^ ln(n)) operations for a matrix of integers of 
fixed size). This algorithm has been implemented using the Maple and 
Giac/Xcas computer algebra systems. 



1 Introduction 



Let's remember that the Jordan normal form of a matrix is 
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where there are 1 or instead of the ?. It corresponds to a full factorization of the 
characteristical polynomial. If the field of coefficients is not algebraically closed, 
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this Jordan form can only be achieved by adding a field extension. The Jordan 
rational normal form is the best diagonal block form that can be achieved over 
the field of coefficients, it corresponds to the factorization of the characteristic 
polynomial in irreductible factors without adding any field extension. 

In this paper, we first present a complex Jordan normal form algorithm. This 
part does not provide an improvement per se, but it gives, in a simpler case, 
a taste of the rational Jordan Normal form algorithm. More precisely we will 
present a similar algorithm that provides a rational normal form maximizing 
the number of Os. This is not a rational Jordan form since the non-diagonal 
block part does not commute with the block-diagonal part, but we show that it 
is fairly easy to convert it to the rational Jordan form. 

This algorithm is not based on the Frobenius form (see e.g. Ozello), and 
assumes that the characteristic polynomial can be fully factorized (see e.g. 
Fortuna-Gianni for rational normal forms corresponding to square-free or other 
partial factorization). It might be combined with rational form algorithm after 
the Frobenius step, but it can be used standalone. It has the same complexity as 
other deterministic algorithms (e.g. Steel), is relatively easy to implement using 
basic matrix operations, and could therefore benefit from parallelism (see also 
Kaltofen et al. on this topic). 

The algorithm of these articles have been implemented in Maple language, 
they work under Maple V.5 or under Xcas 0.5 in Maple compatibility mode. 
They are also natively implemented in Giac/Xcas. Please refer to section 31 to 
download these implementations. 

2 The complex normal Jordan form 

2.1 A simplified case 

Let A be a matrix and -B(A) be (A ■ / — A)^s comatrix. If every eigenvalue is 
simple, we consider one : Aq. Then we can write 

(Ao • / ~ A) • B(Ao) = P(Ao) ■ / = 

The columns of B{Xo) are A eigenvectors for the eigenvalue Aq. To have a base 
of A^s characteristic space for the eigenvalue Ao, we just have to calculate the 
matrix B{Xq) (using Horner's method for example because B{X) is a matrices' 
polynomial) and to reduce the matrix in columns to find one that is not null. 

Our goal is now to find a similar method when we have higher eigenvalues 
multiplicity. 

2.2 Fadeev Algorithm 

First, we need an efficient method to calculate the matrices polynomial B{X). 
Fadeev's algorithm makes it possible to calculate both the characteristic 
polynomial (-P(A) det(A/ — A)) coefficients {pi {i = 0..n)) and the matrices 
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coefficients Bi (i = . . n — 1) of the matrices polynomial giving (A • / — A)'s 
comatrix -B(A). 

{XI - A)B{X) = {XI - A) J2 BkX"" = {J2PkX'')I = P{X)I (1) 

fc<n — 1 fc<n 

By identifying the coefficients of A's powers, we find the recurrence relations : 

Bn-l = Pnl = I, Bk — ABk+i = Pk+ll 

But we still miss a relation between pk and B^, it is given by the : 
Theorem 1 ( Cohen thm ) 

The derivative of the characteristic polynomial P (A), equals the {XI — A) co- 
matrix trace. 

tr{B{X))=P'{X) 

The theorem gives tr(i?fe) = (fc + l)pk+i- If we take the trace in the recurrence 
relations above, we find : 

tr(B„_i) = npn, {k + l)pk+i - tr{ABk+i) = npk+i 

Hence if the field of coefficients is of characteristic (or greater than n) we 
compute Pk+i in function of Bk+i and then Bk : 

tT{ABk+i) n ^ r 

K + 1 — n 
Let's reorder P and B's coefficients : 

P(A) = A"+piA"-^+p2A"-2...+p„ 
B{X) = A"-1/ + A"-2Bi + ...+B„_i 

We have proved that : 

Ai=A, pi = -tr(A), Bi=Ai+piI 

A2=ABi, p2 = --tT{A2), B2=A2+P2l 

Ak = ABk-i, Pk = -■^tr(Afc), Bk = Ak +PkI 

We can now easily program this algorithm to compute the coefficients Bi 
and Pi. The number of operations is 0{n'^) field operations using classical ma- 
trix multiplication, or better 0{n'^~^^) using Strassen-like matrix multiplication 
(for large values of n). For matrices with bounded integers coefficients, the com- 
plexity would be 0{n^ lii('^)) or 0{n"~^'^ ln(n)) since the size of the coefficients 
of Bk is 0(A;ln(A;)). 

Remark 

If tlio field has non-zero characteristic, P(A) should be computed first, e.g. using 
Hessenberg reduction (an O(n^) field operations), then B{X) can be computed 
using Horner division of P(A) by XI— A (an O(n^) field operation using standard 
matrix multiplication) . 
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2.3 Jordan cycles 



Jordan cycles are cycles of vectors associated to an eigenvalue and giving a 
basis of the characteristic space. In a cycle associated to Xq, giving a vector v 
of the cycle, you can find the next one by multiplying {A — Xq ■ I) by v and the 
sum of the sizes of the cycles associated to an eigenvalue is its multiplicity. 

For example, if Aq has multiplicity 5, with one cycle of length 3 and one of 
length 2, the block associated to Aq in the .Jordan basis of the matrix will be : 
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We are looking for vectors giving bases of characteristic spaces associated to 
each eigenvalue of A, and these vectors must form Jordan cycles. 



2.4 Taylor expansion and the characteristic space. 

Let (Ai, m) be the eigenvalues counted with their multipHcities. If the field 
has characteristic 0, we make a Taylor development at the point Ai (cf. equation 
CJ p. El : 

^P(A)/ = {A-XI){B{X,) + B\X,){X-X,) + ... + B^-\X,){X~X,r-') 

= -(A-A,rn(^-^^)"^'^ 

where B'' is the fc-th derivative of B divided by k\. 

If the characteristic of the field of coefficients is not 0, the same expansion 
holds, since the family ((A- Ai)'')*,. is a basis of the vector space of polynomials 
of degree less or equal to rt — 1. In this case (but also in the former case), the 
value of S'"' can be computed using several Horner division of -B(A) by A — Aq. 

As A — XI ^ A — Xil — (A — Ai)/, we have for the m first powers of A — Ai : 

(A-A,/)B(AO = (2) 
{A~XJ)B\X,) = B{X,) (3) 

(4) 

{A-XJ)B^^-\X,) = B^^-\X,) (5) 
(A--A,/)i3"'(A,)-i?"'-^(A,) = -l[{X,-X,r^I (6) 

Theorem 2 The characteristic space associated to Ai is equal to the image of 
S"--i(AO. 

Proof : 

We first show that i3"'~^(Ai)'s image is included in the characteristic space 
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associated to using the fourth equation and the ones before. Let w be a 
vector, V G Jm(S"*~^(Ai)), then 3 m so that v = B"'"^(Aj) • u 



{A-K-ir-v = {A 



= {A~X,-I)-B{X,)-u 

= 

Now we want to prove that every vector v in the characteristic space is also 
in -B"'~^(Ai)'s image. We show it by a recurrence on the smallest integer m 
verifying {A - X,)"'v = 0. 

For m = 0,it's obvious because v — 0. 

Let's suppose the case m is true and prove the case to+ 1. With the equation 
©, we just have to show that : 

w = (A-A,)B"'(A,)i; 

= B''^-\X,)-v-l[{K-X,r^ -v 

is in _B"'~"^(Ai)'s image, because B"'~^{Xi) • u is in B"^~^{Xi)^s image and thus 
if we prove that w is also in, we'll get that ]^(Ai — Aj)"^ • u is in and v is in. 

As B"''{Xi) commutes with A (because it's a polynomial in A) : 

(A - X,rw = B"' (A,) (A - X,r+h^ = 

We can now apply the recurrence hypothesis to w. We now know that w G 

B'''-^{Xi). And so does v. 

2.5 Algorithm 

To find the Jordan cycles, we apply a Gauss reduction on the columns of 
the matrices i?^'^^(Ai) where k < rn. Doing that at the same time for all the 
matrices allow us to keep the relations l(2l to between them after reduction. 

Let's think of the matrices one under another, columns aligned. We reduce 
the matrix B{Xi) and we rewrite the elementary operations on columns done to 
B{Xi) on all the matrices B^{Xi) to keep the relations between them. 

Once the matrix B{Xi) is reduced, if we have k columns not null then we 
know that we already have k rii-long Jordan cycles, taking the k chains of 
vectors associated to the considered column. (In fact /c is or 1 at the first step, 
see the last paragraph in this section). 

If we don't yet have enough vectors to make a base of the characteristic 
space associated to Ai, for each chain of columns of the B''{Xi) corresponding to 
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a column of -B(Ai) that isn't null, wc shift by one matrix down all the columns. 
This process keeps the relations between the matrices. 

Again, We reduce, collect the rij — 1-long Jordan cycles and shift the non-null 
columns. And again as long as we still need vectors to make a basis. 

Remark : 

If there are still columns that are not null after the reduction of -B(Aj), there is 
really only one because one n^-long cycle already gives a base of A^'s character- 
istic space if Aj's multipHcity is Uj. If there are more than one cycle associated 
to Aj, -B(Aj) must be null and we can take care of -B^^^(Ai), etc.. 

2.6 Implementation 

We present here the maple langage implementation. 

2.6.1 Useful functions 

Before implementing this Jordan normal form algorithm, we have imple- 
mented the Fadeev algorithm to calculate the Bi and the characteristic polyno- 
mial's coefficients, then the Horner algorithm to calculate the B{Xi). 

- f adeev(A,Bliste,pliste) takes a matrix A and put the B{X) and char- 
acteristic polynomial's coefficients, in Bliste and in pliste respectively,ordered 
by increasing powers to make it easier to program a polynomial derivation. 
This step requires 0(n'^"'"^) field operations (w = 3 for classical matrix 
multiplication) . 

- evalpolymat(l,a) takes a list of matrices, considered as a polynomial's 
coefficients ordered like before, and a numbera, and gives back the ma- 
trix calculed by the Horner's method. Each evaluation requires 0{n'^) field 
operations (expect n evaluations for a generic matrix with complex coef- 
ficients). 

Then a few utilities : 

- derive_listemat (liste) takes a list as in evalpolymat and gives back 
the derivated list. 

- construction_colonneB(Bliste, pliste) takes what is calulated by fadeev 
and gives back a list of p fists if the characteristic polynomial has p roots. 

In each list, there is first a couple giving an eigenvalue and its multipficity 
and then the matrix of the B^''^ {\i)/k\ for k from to (Aj's multiplicity)—! 
stuck one under another. 

- construct ion (1 ,n) makes the matrix of the eigenvectors using a list of 
eigenvalues and associated cycles, n is the size of the matrix we are study- 
ing, (see the Algorithm part for more details on the fist used by this 
function) . 

The previous section showed that the algorithm requires a reduction in 
columns of the matrix. Maple has a function, called gauss jord, that makes 
reduction but in rows, not columns, so after constructing the column matrix, 
we wfil work with its transposed matrix. To work with it, we needed a few more 
functions : 
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- test_ligne_nulle(B, i) takes a matrix that doen't have to be square, for 
example a matrix with n rows and m columns (we just have to consider 
matrices where n < m) . The function returns 1 if the n first coefficients 
of the i-tli rows are null, if not. 

- decalage_ligne(B, i) takes the partial i-th row (with n coefficients) and 
shifts it right by n. 

- coupe_matrice(B) If the matrix has n rows and m columns ( n < m), 
this function removes the first block nxn. 

2.6.2 The Jordan normal form function 

Splitting the work with all the small functions listed before makes the final 
program quite simple. There are three embedded loops, one loops over all eigen- 
values, it constructs the list of matrices B'^ associated to the eigenvalue, the 
second loop is a while loop that stops when all characteristic vectors for the 
current eigenvalue have been found, the third (inner) loop corresponds to a 
fixed length of the cycles that we are finding. 

The program creates a list of p lists if the matrix has p eigenvalues that are 
all different, each of these p lists contains an eigenvalue and the list of associated 
Jordan cycles. Then with the function described above : construction(l ,n) 
the main program returns the matrix of eigenvectors and the Jordan normal 
form of the matrix A. 

2.6.3 Tests matrices 



A has two eigenvalues : 2 (multiplcity 2) and 1 (multiplicity 1). B has only 
one eigenvalue : 1 (multiplicity 3). Unlike A, the second matrix has two cycles 
associated to only one eigenvalue, it revealed an error in a previous version 
of the program : in the "while" loop, the stop test was inefficient because we 
could collect linearly dependent vectors (because the Maple function "gaussjord" 
making the reduction changes the order of the matrix rows). Hence the test 
function looking if the vector (and the corresponding Jordan cycle we're about 
to collect) is independent of the vectors already collected (by making a matrix 
with all these vectors and searching the rank). 

Once the program showed right for these two examples, it was tested suc- 
cessfully on Jordan matrices constructed with JordanBlock and BlockDiagonal, 
moved to another basis by a random matrix conjugation. 
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2.6.4 Limits of the implementation 

The first version of the program was not really complete because it worked 
only with matrices whose characteristic polynomial, "factors" could factor com- 
pletely (e.g. integer matrices with rational eigenvalues, but not integer matri- 
ces with algebraic eigenvalues). Since "solve" also finds algebraic eigenvalues, a 
"solve"-answer-like to "factors"-answer-like converter was added. Hence this Jor- 
dan normal form program is successfull if and only if "solve" is able to find the 
roots of the characteristic polynomial. 

3 The Jordan rational normal form 

In the previous section, we sometimes had to introduce an algebraic exten- 
sion of the coefficients field (e.g. Q) to be able to compute the characteristic 
polynomial's roots, in this section we will find a basis in the coefficient field 
where the cndomorphism matrix has the best almost diagonal block form, the 
Jordan rational normal form. The diagonal blocks will be companion matrices 
(corresponding to irreducible factors of the characteristical polynomial) , and the 
1 of the complex Jordan normal form will be replaced by identity block matrices. 

We are first going to compute a normal form with as many zeros as possible, 
and from this form, we will compute the Jordan rational form. 

3.1 Pseudo rational Jordan form 
3.1.1 Algorithm 

The method we're going to use is based on an algorithm similar to the one 
used before. Let Q{X) = qo + ■■■ + qd- X'^ be an irreducible factor of the character- 
istic polynomial in the field of coefficients of multiplicity q and degree d of the 
characteristic polynomial P. Note that = 1 since Q divides the characteristic 
polynomial P, hence the euclidean division algorithm of a polynomial by Q does 
not require any coefficient division. 

The characteristic space corresponding to the roots of Q will be replaced 
by a rational characteristic space of dimension d ■ q made of "rational Jordan 
cycles". Recall that : 

{XI -A)- J2 BkX'=P{X)I 

k<n-l 

Since Q{X) ■ I — Q{A) is divisible by A • 7 — ^, there exists a matrix M(A) such 
that : 

{Q{X)I - Q{A))i J2 BuX^)=Q{XyM{X) (7) 

k<n-l 

Now expand B{X) with respect to increasing powers of Q{X) by euclidean divi- 
sion by Q : 

B{X) = J2Ck{X)Q{X)\ deg(Cfe) < q 

k 
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Replacing in and observing that the matrix coefficients of order less than d 
vanish, we get : 

QiA) • Co = 0, Cfe = Q{A) ■ Ck+i 

This is similar to the case where the eigenvalue is rational, we get a chain of 
polynomial matrices that are images of the preceding one by Q{A) : 

Cq-l — ^ Cq-2--- Co ^ 

We will find the rational Jordan cycles by constructing Jordan cycles of Q{A). 
Note that if we find a Jordan cycle of length k for Q{A) we can construct d — 1 
other Jordan cycles by multiplying the cycle by A^ for i = l..d — 1. 
All these vectors are independent, indeed if 

kjA'Q{Ayv = 0, Q{A)''v = 0, QiA)''-^v + 

by multiplying by Q{A^^~^ we get : 

\,u^^A')Q{Af-\ = ^ 

i 

hence \i^k-i — for all is since Q{A)^^^v ^ and Q{A) is irreducible. Multi- 
plying further by Q{A)^^'^ , identity, it follows that all A^.j are zero. 

Once we have collected these kd vectors, we search for another cycle in the 
vectors of the Cj matrices that are linearly independant to all A^Q{A)^~^v 
starting from Co and increasing j. If we find a new end cycle vector Q{A)^ '^w 
such that Q{A)^ w — Q and Q{A)^ is independent of the preceding end-cycle 
vectors, then we can form k' d vectors A^Q{Ayw. We will show that these vectors 
are independent of the A^Q(Ayv since (3(A) = go + •■ + 9d • '^d is irreducible. 
Indeed if we had a relation like 

KoA'Q{Ayv + ^l,^JA'Q{Ayw = o, 

If j > k' then Xij = by multiplication by Q{Ay for decreasing j > k' . Now 
we multiply by Q{A)'' and we get two polynomials P and R of degree less 
than degree((5) such that : 

P{A)Q{A)''-\ + R{A)Q{Af~^w = 

Since Q is irreducible, it is prime with i? if i? 7^ 0. Hence if i? 7^ 0, by applying 
Bezout's theorem, we could invert R modulo Q and express w as a linear com- 
bination of A'Q{A)''~'^v. Therefore i? = and P = and Hi,k'~i = K,k'-i = 0. 
Let (vk-i) {vk-2) ■■■ (vq) (0) be a cycle of Q{A), we have : 

{vk-i,Avk-u-,A''-^Vk-i) ... ^ {vo,Avo,...,A'^^'^vo) -> (0,...,0) 

where the arrow means "image by Q{A)". 
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Let's write the matrix A in the base fo, Avq, .., A'^ ^vq, ■■, Vk-i, ■■, A''' ^Vk-i 
we find an "almost Jordan rational blockl", its size is k ■ d : 
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Indeed vq image hy A is A ■ vo the second vector basis, etc. to A''' ^ ■ vo whose 
image by ^4 is : 

A'^-Vo = (QiA) -qo-q,.A-...- Qd-i ■ A''-^) ■ vo 

Since Q{A) ■ vq = {vq ends a Jordan cycle of Q(A)), we get the first block of 
the matrix in the new basis. 

For the second block, we get the first d—1 columns in a similar way. For the 
last one : 



A'^-vi = {Q{A) -qo-qi-A 



Qd-l ■ A""-') ■ VI 



Since Q{A) ■ vi ~ vq, we get the above matrix part. By applying the same 
method to the rest of the cycle we get the matrix. 



3.1.2 Complexity 

Each euclidcan division requires 0{n^d) field operations (d is the degree of 
the irreducible factor) . There are q euclidean divisions of a polynomial of degree 
less than n with n, n matrices coefficients by a polynomial of degree d, hence 
computing the Cj requires 0{n'^dq) operations, adding for all irreducible factors, 
we get a complexity of 0{n'^) for the division part. 

Let rid,...,rqd be the number of Jordan cycles of Q{A) of length q, 1. 
We have : 

riq + r2{q - 1) + ... + = g 

The first step of the reduction part requires reducing a n, nq matrix of rank nd. 
Then we will reduce a rid+n, n(g — 1) matrix of rank {ri+r2)d such that the rid 
first rows are already reduced and independant (hence r2d new independent rows 
in the n last rows remain to be extracted), etc., then a {ri + ... + ri)d+n,n{q — i) 
matrix of rank (ri + ... + ri_|_i)d with first (ri + ... + ri)d independent reduced 
rows and rt+id new independent rows in the n last rows to extract, etc. We 
will have to make nrid row operations on the i-th matrix. Hence we will make 
0{nridn{q — ij) operations on the i-th matrix. Adding all reduction steps, we 
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will make Oijn^dq) field operations for each irreducible factor, hence 0{n^) field 
operations for all irreducible factors. 

The complexity of the whole pseudo-rational form is therefore 0{n^) field 
operations and is dominated by the Cj computation (since B can be computed 
in 0{n'^'^^) field operations). 

3.1.3 Example 
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The characteristic polynomial of A\s{x — 2)^(a;^ — 2)^. For A = 2 there are 2 
eigenvectors : 
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For a;^ — 2 of multipHcity 2, we find a cycle of length 2 for Q{A) = A? — 2- 1 : 

(0, 0, 0, -1, -1, -1) ^ (1, 0, 0, -1, -1, -1) ^ (0, 0, 0, 0, 0, 0) 
After multipHcation by ^4, we get : 

((0, 0,0, -1,-1,-1), (1,4, 1,4, 0,-3)) ^((1,0,0, -1,-1,-1), (2,4, 2,4,0, -2)) -^0 
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The matrix P is therefore : 
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To obtain the rational normal form, we must replace the block 
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3.2 From pseudo-rational to rational Jordan form 

The pseudo rational form has unfortunately not the commutation property, 
the block diagonal part does not commute with the remainder, hence we will 
compute the rational Jordan form from the pseudo rational form. 

We now assume that we are in a basis where the endomorphism is in pseudo 

/... 1\ 

rational form, and we want to compute a new basis so that the I ... I 

blocks are replaced by identity matrices. Let's assume that we have made the 
first j blocks (each of size d) indexed from to j — 1 corresponding to the family 
of vectors (t;o,0) •••) I'Ojd-i) •••) We want to find a vector vjfi to begin 
the next block. The Vj^i will be defined in function of vj^i-i using the relation 
Avj^l-i = Vj^i + Vj-i^i-i. Hence vjfi must satisfy : 

Avj^d-i = -QoVjfi - ... - qd-iVj,d-i + Vj-i^d-i (8) 
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Applying the previous recurrence relations, we determine Q{A)vjfi with respect 
to Vjij (with j' < j, I < d). Since Q(A) is a shift of d indices to the left, we 
will let Vjfi be the shift of d indices of Q{A)vjfi to the right (if we stay in the 
original basis, "inverting" Q{A) can be done using the pseudo-rational basis). 

More precisely, let's compute vjj in terms of the Vj^ and Vj\i' {j' < j). We 
denote the binomial coefficients by (5„) (they can be computed efficiently using 
Pascal's triangle rule). A straightforward recurrence gives : 

inf(ij) 

Vjj = A^Vj^o - ('^) Vj_rn,l-m (9) 

m—1 

Replacing in we get : 

inf(£ij) d-l mf(;j) 

m—1 1=0 m—1 

eventually : 

d inf(ij) 

m) '^j — ra.l—ra (10) 

1=1 m=l 

Application to the example : 

We stay in the original basis for the coordinates. Here uo,o — (4, 24, 12, 32, 8, —4) 
and -wo.i — Avjfi. A preimage by Q{A) is given by wi^ = (0, 4, —4, 8, 4, —4) and 
wi^i = Awifi. Applying ljl()|l . and 51 = 0, (72 = 1 we must satisfy : 

2 inf(i,i) 

Q{A)vifi = ^ g; Y (m) = 2^^0,1 

1=1 m=l 

hence : 

wi,o = 2^(0,4,-4,8,4,-4) = (-8,-32,0,-48,-16,16) 
vi,i = At;i,o-wo,o - (4,40,-4,64,24,-20) 

We have indeed Avi^i = 2vifi + wo,i. 

3.3 Maple implementation 

In the first part, we were working with matrices polynomials and not poly- 
nomial matrices, so the first thing to do was to create a traduction function 
(which takes the list of the matrices that are S's coefficients, B given by Fadeev 
algorithm) to make the euclidean divisions on B{X) coefficient by coefficient 
(nouvelle_ecriture function, arguments are B the polynomial Q we want to 
divide by, and Q's multiplicity). Then we collect the cycles of Q{A) as in the 
complex Jordan form case, by gluing the d matrices vertically and transposing 
the result for Gauss- Jordan reductions. The main changes are that we generate 
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cycles of A by multiplication by I, A, A"^'^ (f abriq_cycles function) and 
we must take care that a new end-cycle vector must be independent not only of 
a previous end-cycle vector Vi but also of its images {A ■ Vi, A'^~^ ■ Vi}. 
The structure of the main rational Jordan form function Jordan2 is : 

- A call to demarrage that will return a list of [[irreducible polynomial, 
multiplicity], [cycles]]. 

- For first order irreducible polynomials, the functions of the complex nor- 
mal form are called 

- For each irreducible polynomial, conversion from pseudo-rational Jordan 
form to rational Jordan form 

- a call to construction_special to build the passage matrix. 

4 "User guide" 

The Giac/Xcas free computer algebra system is available at : 
www-f ourier . uj f -grenoble . f r/~parisse/giac . html 

The functions Jordan and rat_ Jordan implement the Jordan normal form and 
the rational Jordan normal form. 

The maple implementation of this algorithm is available at : 
www-f ourier .ujf -grenoble . fr/~parisse/ Jordan. map 

Once the Maple session is opened, run the command readC'Jordan.map") • 
Then three programs are available : 

- TER_ Jordan takes a matrix A and returns the matrix of eigenvectors and 
the Jordan normal form of A. 

- final takes a matrix A and returns the matrix of eigenvectors and the 
pseudo-rational form, calculated with a hybrid method combining the two 
programs above. 

- Jordan2 takes the matrix A and returns the rational form. 

Note that in the current version, there is a small inconsistency, since for the 
rational roots of the characteristical polynomial, the Jordan 1 are not on the 
same side of the diagonal than the Jordan identity blocs for irreducible factors 
of degree larger than 1. 

This Maple implementation can also be run under Xcas, but it is of course 
much faster to call the native Xcas functions. 
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